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Abstract 

We present a new class of exact solutions for the so-called Lapla- 
cian Growth Equation describing the zero-surface-tension limit of a 
variety of 2D pattern formation problems. Contrary to the common 
belief, we prove that these solutions are free of finite-time singulari- 
ties (cusps) for quite general initial conditions and may well describe 
real fingering instabilities. At long times the interface consists of N 
separated moving Saffman- Taylor fingers, with "stagnation points" in 
between, in agreement with numerous observations. This evolution 
resembles the N-soliton solution of classical integrable PDE's. 
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The problem of pattern formation is one of the most rapidly developing 
branches of nonlinear science today. [1] Of special interest is the study of 
the front dynamics between two phases (interface) that arises in a variety 
of nonequilibrium physical systems. If, as it usually happens, the motion 
of the interface is slow in comparison with the processes that take place in 
the bulk of both phases (such as heat-transfer, diffusion, etc.), the scalar 
field governing the evolution of the interface is a harmonic function. It looks 
natural then, to call the whole process Laplacian growth. Depending on the 
system, this harmonic scalar field is a temperature (in the freezing of a liquid 
or Stefan problem), a concentration (in solidification from a supersaturated 
solution), an electrostatic potential (in electrodeposition), a pressure (in flows 
through porous media), a probability (in diffusion-limited aggregation), etc. 

We present in this paper a new class of solutions of the 2D-Laplacian 
growth problem in the limit of zero surface tension. These solutions are 
quite general because no symmetries of the moving interface are assumed. 
Most remarkably, they do not develop finite time singularities but, contrary 
to the common belief, remain smooth for all finite times. Thus, they may 
describe real fingering instabilities when surface tension is very small. In 
the long-time limit they give rise to N separated fingers, each of which (for 
enough separation) describes the Saffman- Taylor finger [2] in channel ge- 
ometry, and whose evolution closely resembles the iV-soliton formation in 
nonlinear integrable PDE's. 

In the absence of surface tension, whose effect is to stabilize the short- 
wave perturbations of the interface, the problem of 2D Laplacian growth is 
described as follows: 




(1) 
(2) 
(3) 
(4) 



d„u\j: = 1, 



Vn = -dnU\r{t)- 
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Here u{x, y; t) is the scalar field mentioned above, T{t) is the moving interface, 
E is a fixed external boundary, (9„ is the component of the gradient normal 
to the interface (i.e. the normal derivative), and Vn is the normal component 
of the velocity of the front. 

We consider in this paper an infinitely long interface. We introduce then 
a time- dependent conformal map / from the lower half of a "mathematical" 
plane, C, = ^ + ir], to the domain of the physical plane, z = x + iy, where 
the Laplace equation (1) is defined, ( ^ z. We also require that f{t, C) ~ C 
for C — > ^ — ioo. With this definition, the function z = f{t,^) describes the 
moving interface. 

Using this conformal map and taking into account the boundary condi- 
tions of the problem we find: 

- (^) 

and 

where the overbar means complex conjugate, the subscripts t and ^ indicate 
partial derivatives with respect to t and ^, respectively, di is the component 
of the gradient tangent to the interface, ip is a harmonic function of x and 
y, conjugate to u, that satisfies ip = due to the boundary conditions (2) 
and (3). If we equate Eqs. (5) and (6), in accordance with Eq.(4), we finally 
obtain: 

Mftfd = 1 (7) 

As in Ref. 3, we will refer to Eq. (7) as the Laplacian Growth Equation 
(LGE), because the scalar field determining the growth obeys the Laplace 
equation (1). The LGE was first derived, to our knowledge, in 1944 in- 
dependently by Polubarinova- Kochina[4] and Gahn[5]. This equation has 
remarkable properties, unexpected for nonlinear PDE's, such as an infinite 
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set of conserved quantities [6], a pole decomposition [7], an exact solution in 
the class of finite polynomials [8] and other exact solutions (though for partic- 
ularly symmetrical initial shapes) [9]. Also the Saffman- Taylor finger [2] is a 
particular travelling wave- solution of this equation. All these properties (ex- 
cept the latter one) are nontrivial and nonperturbative due to the nonlinear 
nature of the LGE. 

Unfortunately, despite these remarkable properties, practically all known 
solutions of the LGE show finite-time singularities via the formation of cusps 
[7], [8], [10]. Therefore, all these nonperturbative results are helpless to shed 
light on the physics and geometry of the system in the long-term limit. (Al- 
though a few exact results have been presented that have no finite-time cusps 
[9], they correspond to cases with very restricted symmetries of the moving 
interface.) As a result it has been generally assumed that these finite time 
singularities arc an essential feature of LGE solutions [11] and that, in this 
sense, the physics represented by the LGE is incomplete. Thus, the natural 
attitude was to include surface tension in the theory to stabilize the moving 
interface and get rid of the finite-time singularities. [11] 

The main result of this paper is to show that, contrary to this widespread 
view, the LGE (i.e. a zero-surface tension limit of the 2D Laplacian pat- 
tern formation) admits quite a broad class of exact time-dependent solutions 
which remain smooth for an infinite time. These solutions arc of 'W-finger 
type" because they lead to the formation of N well-developed and separated 
fingers in the long-time asymptotics (see Fig.l). Thus, the problem of finite- 
time cusps can be solved within the framework of the LGE. (However, there 
are problems in which the inclusion of surface tension is indeed unavoidable, 
e.g. the selection mechanism in dendritic growth). [1] 

To introduce this new class of solutions, we start from the statement that 
any function f{t, () whose derivative, has an arbitrary distribution of 
moving poles, Ck{t) = + ''■Vk, and roots, Zk{t), in the upper-half plane, Im 
C > 0, and no other singularities, is a solution of the LGE. This is easy to 
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verify by substitution of any such function / into Eq. (7). Let us consider 
the arbitrary case when has A'" + 1 simple poles. Since we assume that 

with all poles and roots initially lying above the real axis, we have 

Af+l 

f = A{t) + C-iT. aUog(C-a) (9) 

k=l 

where A{t) is an unknown function of time, and the factor —i is taken just 
for convenience. 

By substitution of the last equation for C, real, into the LGE (Eq. (7)), 
we find that 

• i) A{t) = -it ; 

• ii) ctfc G C does not depend on time ; 

• in) the time-dependence of the poles Cfe is governed by the following 
constants of motion 

Af+l 

(^k = f{Ck)^-it + Ck-iY. «^log(a-0), l<k<N + l. (10) 

t=l 

Therefore, the motion of the interface is determined by Eqs. (9)- (10) for real 
C^. This solution, for N = 1 corresponds to the development of an isolated 
finger, similar to the one found by Saffman for channel geometry. [12] 

The break of analyticity of the interface (a cusp) occurs when at least one 
of the moving poles, Cfc(^)i or zeros, Zk{t), of crosses the real axis, = 0, 
of the mathematical plane, (. If all Cfc's and Z^s remain on the upper-half 
plane during the whole evolution, then the moving interface remains smooth 
(analytic) for an infinite time. To obtain sufficient conditions under which 
this is true for / given by Eqs. (9)-(10), we note the following: 
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(i) In order for the solution to exist as t ^ oo and satisfy r]k > for all 
finite times, all a^'s must have positive real part. This is the only way that 
the divergent term —it in the R.H.S. of Eq. (10) can be compensated. The 
term that compensates it is 

- iRe{ak) log(Cfe - Cfe) = -iRe{ak) \og{-2ir)k) (11) 

and implies that r^^ — as t — cxo. 

(ii) An isolated singularity, (k, can never reach the real axis at a finite 
time. If it did, the term —iaklog{(k — Ck) in Eq. (10) would diverge and 
could not be compensated by any other. On the other hand, if all Read's 
are positive, then groups of M < + 1 singularities could not reach the real 
axis simultaneously at a finite time for exactly the same reason. Therefore, 
we can conclude that, if Rea^. > and r]k{t = 0) > for 1 < A; < + 1, then 
r)k{t) > for all finite times t and 77^ as t — cxd. 

(iii) We assume first that all a^s arc real and positive. We calculate the 
derivative of / with respect to ^ = ^ + ir;. After a little algebra, we write its 
real part as 

Re/c = l+E^:^-^E^^ (12) 

k=l Is Sk\ IS Sk\ 

Since all a^s and all rik{t = 0)'s are positive, then by the result in (ii), also 
rjk > for all finite times. Thus, we see from Eq. (12), that Re/^ equals zero 
only if 1] is strictly positive. This means that the zeros, Z^, of lie always 
on the upper-half of the mathematical plane and never cross the real axis. 

Consequently, if at time t = 0, all 's are real and positive and all 
rik's are positive, then the interface represented by Eq.(9) remains smooth 
throughout its evolution. 

For the general case of complex a^s (with positive real part) we are 
currently unable to prove neither the existence nor the absence of finite-time 
blow-ups. (For A^ = we can easily show that the solution (9) remains 
smooth for all time). In spite of the lack of a rigorous proof, we believe that 
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also in this case, there are no finite-time cusps for a broad range of initial 
conditions. This conjecture is supported by a series of numerical experiments 
that we will discuss in a forthcoming paper. [13] 

Eq. (9) is not the only solution of the LGE that is characterized by the 
motion of simple poles. For example, if in Eq. (9) we replace log{( — (k) 
by log(e*'' — 6*''*=), and introduce a new parameter A, we find a 27r-periodic 
solution, relevant for channel geometry, of the form: 

t ^ 

f = -i + - iJ2 ak log{e"^ - e'^'^^ - tt < Re C < tt. (13) 

^ k=i 

This sohition must satisfy J2k ctfc = 1 — A so that = 1 for C ^ —too. The 
parameter A is the fraction of width of channel occupied by the fingers. This 
solution has the same properties as Eq. (9). In particular, there exist 
constants of motion defined by /4 = f{(k) and cusps are absent if all a^s 
are real and positive and all rjkS are positive, as it follows from the equation: 

N 

/C - A + ^ e»(l-0|2 (^ - e-'''e'^<-^^^) ^ if ImC < 0. (14) 

k=l I I 

Eq. (13), but with e~*^ instead of e"", was proposed as a solution of the LGE 
in [11], where cusps were found via numerical simulations. Taking advantage 
of the corresponding constants of motion, we can easily show the necessity 
of these cusps if Re{ak) > (as is the case in [11], where — I for all k). 

The constants have a clear geometrical meaning in the physical plane 
for both solutions (9) and (13): |7rQ;fc| is related to the width of the gap 
between adjacent moving fingers (in the case of enough separation), while 
arg(Q;fc) is related to the angle that this gap forms with the horizontal. [13] 
We show this property in Fig. 1, where we have plotted the interface Y = 
Im fit, x) vs. X = Ref{t,x), for the solution (9) with two singularities 
{N = 1), at a particular time. The real parts of both a's are positive, but 
while CKi is purely real (ai = 0.8), a2 has a nonvanishing imaginary part (02 = 
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0.8 + iO.l). We have drawn on top of each gap a dashed hne of length {nakl 
and slope lm{ak) /Re{ak) that highlights the meaning of at- We conclude 
then that, if all a^'s are real and positive, the interface is a single- valued 
function Y{X). Among the experiments on two-dimensional viscous fingering 
in a channel, one can find different degrees of bending and ramification of 
the moving fingers. [14] It follows from these experiments that non-single- 
valued interfaces generally appear, although with small bending. Therefore, 
complex a's (possibly with small imaginary parts) would be necessary to 
describe them using Eq. (9). As we mentioned before, we expect that even 
in this more general setting, Eq. (9) still provides a meaningful (non- singular) 
solution, especially, if the imaginary parts of a's are small enough. 

The constants f3k also have a clear geometrical meaning in the physical 
plane: the points (Re/?fc! Iiii/?^ + Re log 2) are the coordinates of the tips 
of the gaps between fingers. We show this property in Fig. 2, where we have 
plotted the interface obtained from Eq. (9) at different times for a particular 
choice of + 1 = 7 singularities and real and positive a^s. For the sake 
of generality we have deliberately taken the initial condition without any 
particular symmetry. We have indicated in this figure the location of the tips 
with asterisks. As one can see, they are "stagnation points" of the interface. 
This kind of stagnation points have been observed in numerous experiments 
[14] and numerical simulations [15]. Because all real experiments deal with 
non-zero surface tension, the experimental evidence of stagnation points in 
viscous fingering supports our conjecture that a small surface tension does 
not destroy the behavior of the zero-surface-tension solutions. It is interesting 
to note that we can also identify stagnation points occurring in numerical 
[16] and physical [17] experiments of diffusion-limited aggregation (DLA), 
(a feature usually explained as a "screening" effect). This is in accordance 
with the view (which we believe) that DLA fractal growth is described by 
the LGE in the continuous limit. [18] 

Finally, we show in Fig. 3 the motion of the singularities in the mathe- 
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matical plane that correspond to the solution depicted in Fig. 2. Comparing 
both figures we can see that, while the singularities tend to the real axis, 
the interface takes the shape of a set of more and more well-developed and 
separated N fingers. The evolution of this interface resembles the A^-soliton 
solution of classical exactly solvable PDE's, such as the Korteweg-de Vries 
or nonlinear Schrodinger equations, where in the long-time asymptotics one 
can also have separated solitons, each of which described by the single 
soliton solution of the corresponding PDE. An evident difference is that fin- 
gers, unlike "classical" solitons, always have a non-zero velocity component 
normal to the interface. The connection between the N-soliton and the N- 
finger solutions is deeper than a superficial resemblance, as we will show in 
a forthcoming paper. 

The separation of the singularities in the long-term asymptotics and the 
similarities with the A^-soliton solutions suggests that there might be some 
transformation into action-angle variables in which this separation arises nat- 
urally. We also suspect that some Hamiltonian structure (at least formal) 
might underlie the whole evolution. We arc presently working in this direc- 
tion. The nontrivial extension of this work to the case of a closed (finite) 
interface is also of interest. In this case, our proof regarding the absence 
of finite-time singularities does not hold, however, the A^-finger like solution 
with its associated constants of motion is still valid. The observation in 
physical and numerical experiments of very similar behavior to the infinite 
interface, including the existence of "stagnation points" [19], makes it seem 
likely that the class of solutions studied in this paper may also shed light on 
the problem of radial geometry. 

We thank M. Ancona and R. Camassa for their useful comments. This 
work is supported by the US Department of Energy at Los Alamos National 
Laboratory. 
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Figure Captions 



Figure 1. We plot the interface, Y =Im /(t, as a function of X =Re 
f{t,^), at a particular time, for / defined in Eq. (9) with A — —it and two 
singularities, ai — 0.8; 0:2 = 0.8 + iO.l. The dashed lines that we have drawn 
on the gaps between fingers have length |7rQ;fc| and slope lm{ak)/Re{ak). 
Figure 2. We plot the interfaces, y as a function of X, at times t — 0, 
f = 9, i = 24, and t = 30, for / given by Eq. (9) with A{t) = -it, ( 
real and N + 1 — 7. The evolution of the seven singularities was obtained 
using a fourth-order Runge-Kutta integrator to solve the set of evolution 
equations that are obtained by differentiating with respect to time Eq. (9) 
with A{t) = —it. The paramaters of the simulations are ai — 0.80, Pi — 
6.00 - i37.n,a2 = 2.00 P2 = 20.44 - i33.98, as = 1.00, P3 = 35.61 - i34.49, 
q;4 = 0.50, /34 = 43.15 - i36.31, 0:5 = 1.50, 05 = 54.58 - i32.78, ag = 0.35, 
Pe = 64.65 -i40.88, ^7 = 1.80, and P7 = 72.90 -i37.40. We see that, from an 
initially bell-shaped form, the interface evolves into a situation characterized 
by = 6 well separated fingers. We also observe the existence of "stagnation 
points" of the interface (indicated with asterisks) . 

Figure 3. We plot the trajectories of the moving singularities, Cfc(^)) on the 
mathematical plane, that correspond to the physical process dcspicted in Fig. 
2. The asterisks indicate, from top to bottom, the points on the trajectory 
at times t = 0, t = 9, t = 15, and t = 30. Wc see that all the singularities go 
asymptotically in time to different values and that all 77 — > as t — > 00. 
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